Fluctuations and vortex pattern ordering 
in fully frustrated XY model with honeycomb lattice 



S. E. Korshunov 1,2 and B. Doucot 2 

1 L. D. Landau Institute for Theoretical Physics, Kospgina 2, Moscow 119334, Russia 
2 Laboratoire de Physique Theorique et Hautes Energies, CNRS UMR 7589, 
Universite Paris VI and VII, 4 place Jussieu, 75252 Paris Cedex 05, France 
(Dated: February 2, 2008) 

The accidental degeneracy of various ground states of a fully frustrated XY model with a hon- 
eycomb lattice is shown to survive even when the free energy of the harmonic fluctuations is taken 
into account. The reason for that consists in the existence of a hidden gauge symmetry between 
the Hamiltonians describing the harmonic fluctuations in all these ground states. A particular vor- 
tex pattern is selected only when anharmonic fluctuations are taken into account. However, the 
observation of vortex ordering requires relatively large system size i>I c > 10 5 . 
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A fully frustrated XY model can be denned by the 
Hamiltonian 

H = - J y^cosjcpj -(fj- Ajj) , (1) 

m 

where J > is the coupling constant, the fluctuating 
variables <p+ are defined on the sites i of some regular 
two-dimensional lattice, and the summation is performed 
over the pairs of nearest neighbors (ij) on this lattice. 
The non-fluctuating (quenched) variables Ay = — Aji 
defined on lattice bonds have to satisfy the constraint 
^ Aij = 7T (mod 2tt) (where the summation is performed 
over the perimeter of a plaquette) on all plaquettes of the 
lattice. 

For two decades such models (on various lattices) have 
been extensively studied Q in relation with experiments 
on Josephson junction arrays Q, in which ipi can be as- 
sociated with the phase of the superconducting order pa- 
rameter on the i-th superconducting grain, and is 
related to the vector potential of a perpendicular mag- 
netic field, whose magnitude corresponds to a half-integer 
number of superconducting flux quanta per lattice pla- 
quette. Planar magnets with odd number of antiferro- 
magnetic bonds per plaquette are also described by 
fully frustrated XY models. Recently, the active inter- 
est in fully frustrated Josephson arrays has been related 
to their possible application for creation of topologically 
protected quantum bits 

The ground states of the fully frustrated XY models 
are characterized by the combination of the continuous 
U(l) degeneracy (related with the possibility of the si- 
multaneous rotation of all phases) and discrete degener- 
acy related with the distribution of positive and negative 
half-vortices between the lattice plaquettes. Since vor- 
tices of the same sign repel each other, the energy is 
minimized when the vorticities of the neighboring pla- 
quettes are of the opposite sign. In the case of a square 
lattice this requirement is fulfilled for all pairs of ncigbor- 
ing plaquettes when the vortices of different signs form 



a regular checkerboard pattern Analogous pattern, 
in which the vorticities of the neighboring plaquettes are 
always of the opposite sign, can be constructed in the 
case of a triangular lattice ja, Q ■ 

In the case of a honeycomb lattice it is impossible to 
construct a configuration in which the vorticities are of 
the opposite sign for all pairs of neighboring plaquettes. 
As a consequence, the discrete degeneracy of the ground 
state turns out to be much more developed 0,0] , and can 
be described in terms of formation of zero-energy domain 
walls in parallel to each other ^(Ji the residual entropy 
of the system being proportional to its linear size. Quite 
remarkably, the comparison of the free energies of weak 
fluctuations in two different periodic ground states shows 
that in harmonic approximation they are exactly equal 
to each other Jl^}, although the spectra of fluctuations 
A(k) in these states are essentially different: for example, 
for small momenta k they are characterized by different 
values of (anisotropic) helicity moduli. 

In the present work we demonstrate that this ab- 
sence of degeneracy lifting is not a simple coincidence, 
but a consequence of a hidden gauge symmetry between 
the Hamiltonians of harmonic fluctuations in different 
ground states, and extends itself to all states formed 
by some sequence of parallel zero-energy domain walls. 
Therefore, the inclusion of harmonic fluctuations will not 
lead to removal of the accidental degeneracy even if in- 
stead of considering thermodynamic fluctuations one cal- 
culates the free energy of quantum fluctuations at arbi- 
trary (or zero) temperature. 

This gauge symmetry is broken when anharmonicities 
are taken into account. We have compared the lead- 
ing anharmonic contributions to the free energies corre- 
sponding to different ground states and have found which 
of these states is selected at low temperatures. How- 
ever, the difference between their free energies turned 
out to contain extremely small numerical coefficient. 
As a consequence, the observation of vortex ordering 
is possible only in relatively large systems (whose size 
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L>i c > 10 5 ). 

We believe that the discovery of these new features 
(some of which are completely unexpected) makes our 
results of interest not only in relation with Josephson 
junction array physics, but also in the more general con- 
text of two-dimensional statistical mechanics. 

Fig. la shows a structure of the simplest ground state 
of the fully frustrated XY model with a honeycomb lat- 
tice. Each arrow corresponds to 0y = (pj — <pi — Aij = 
7r/4, whereas on the bonds without arrows 9ij = 0. In 
this state the plaquettes with positive and negative vor- 
ticities form straight stripes, so in the following we shall 
call it a striped state. 

A striped state allows for formation of domain walls 
(separating two different realizations of such state) which 
cost no energy |ldl | . An example of such zero-energy do- 
main wall is shown in Fig. lb. An arbitrary number 
of zero-energy domain walls separated by arbitrary dis- 
tances can be introduced into the system in parallel to 
each other [Io|. 

If such domain walls are created at every possible po- 
sition, another periodic ground state is obtained, which 
is shown in Fig. lc. In accordance with the shape of the 






FIG. 1: Three different ground states of a fully frustrated XY 
model with a honeycomb lattice: (a) striped state, (b) zero- 
energy domain wall and (c) zig-zag state. Directed arrows 
(or their absence) correspond to Oij = 7r/4 (9ij = 0). Lattice 
plaquettes with positive vorticities are shaded. 



lines formed by the plaquettes with positive and negative 
vorticities in this state we shall call it a zig-zag state. 
Alternatively, one can describe all other ground states as 
obtained by the introduction of zero-energy domain walls 
on the background of a zig-zag state. 

If all sites of a honeycomb lattice are numbered by pairs 
of integers (n,m) as shown in Fig. 2, the Hamiltonian 
describing the harmonic fluctuations in the striped state 
of Fig. la can be written as: 



H 



(2) 



n m=n(mod2) 



[Ml 



+ JliUn.m — Vn+1 f+JiiUr, 



>-i) 2 + (2) 
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^n— 1,771 ) 



where J\ = J2 = Jcos(7r/4), J3 = J, whereas u n ^ m and 
v n ,m are the deviations of the variables tp from their equi- 
librium values on two triangular sublattices forming a 
honeycomb lattice. 

If one assumes the presence of periodic boundary con- 
ditions in the horizontal direction and open boundary 
conditions in the perpendicular (vertical) direction, the 
introduction of plane waves with respect to the variable 
n allows to rewrite Eq. (J2J as: 



dq 



Y,{Jsl\u m (q)\ 2 + \v m ( q )f 



2tt 

m 

Ji[u m+ i(q)v* m {q) 



(3) 



[K{q)u m {q)v* m {q) + c.c.]} 



where J$ = J\ + J2 + J3 and K(q) = J2 cxp(— iq) 
J^eKpiiq). A trivial gauge transformation: 



u m +i(q) 
v m (q) 



exp[ia(q)m] x 



u m+ i(q) 
v m (q) 



(4) 



where a(q) = arg[K(q)], allows then to replace Eq. © 
by the analogous expression with Ko(q) = \K (q)\ substi- 
tuted for K{q). 

Note that in the chosen mixed representation the only 
modification of the Hamiltonian which appears when an 
arbitrary sequence of horizontal domain walls is intro- 
duced consists in replacement of K (q) by K*(q) for some 
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FIG. 2: The numbering of a honeycomb lattice sites by a pair 
of integers (n,m) used when writing Eq. @. 
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values of m. It is rather evident that the gauge trans- 
formation analogous to Eq. (in which a(q)m should 
be replaced by a(q) J2 m '<m s ™'' wnere variable s m = ±1 
describes the choice between the two options existing for 
the continuation of a ground state at each m) allows to 
transform any such Hamiltonian to the same form (with 
K{q) replaced everywhere by Ko(q)). 

This means that for the boundary conditions described 
above, the whole set of eigenvalues will be exactly the 
same for all Hamiltonians obtained by the introduction 
of an arbitrary sequence of horizontal domain walls (even 
for a finite sample). Accordingly, the free energy of the 
harmonic fluctuations will be exactly the same, and can- 
not be the source for the selection of a particular ground 
state. 

Clearly the free energy of harmonic fluctuations also 
remains degenerate when one considers a quantum gener- 
alization of the same model with the diagonal mass term, 
which in terms of a Josephson junction array corresponds 
to taking into account the self-capacitance of each super- 
conducting island . In order to include into considera- 
tion the mutual capacitances of neighboring islands (that 
is the capacitances of the junctions), one has to apply the 
same approach (construction of the gauge transformation 
which makes all the coefficients real) not to the harmonic 
part of the Hamiltonian, but to the frequency dependent 
Fourier transform of the harmonic contribution to Eu- 
clidean Lagrangian, which also turns out to be possible. 
That means that in the quantum version of the model the 
accidental degeneracy survives (at the harmonic level) for 
arbitrary relation between the self-capacitance of an is- 
land and the capacitance of a junction. 

If one assumes now periodic boundary conditions in 
the vertical direction (instead of open boundaries), the 
degeneracy of the free energy associated with harmonic 
fluctuations (quantum or thermodynamic) will be mani- 
fested only in the thermodynamic limit. 

The accidental degeneracy of different ground states 
is removed when anharmonicities are taken into account. 
The leading contribution to the free energy induced by 
anharmonic fluctuations, F& n h, is given by the sum of 
two terms, which in the classical limit can be written 
as = -([#( 3 )] 2 )/2T and = (ffW), where 
and ij' 4 ' are, respectively, the third- and the fourth-order 
corrections to the harmonic part of the Hamiltonian. 

Since each term in F^ depends only on local phase 
difference on a particular bond, it can be proven with the 
help of the hidden gauge symmetry discussed above that 
in the considered system F^ is the same for all ground 
states. However, this property does not extend itself to 
F^ 3 \ which depends also on more distant correlations. 

Comparison of the expressions for F^ in the two dif- 
ferent periodic ground states shows that the main con- 
tribution to SF = _ — Fjjfa (normalized per single 



hexagon) can be written as 

6F=-6(G 3 S + G S Z )J*/T , (5) 
where J* = Jsin(7r/4)/6, whereas 

G s = ({u n ,m ~ Vn+l,m)(u n +l,m-l — W n +2,m— l)) (6) 

and 

G z = ((w„,m - l'n+l,m)K+l,m-l ~ Un.m-l)) (7) 

are two correlation functions (for the bonds with 

- calculated with the help of the harmonic Hamilto- 
nian in striped and zig-zag states respectively. 

Numerical calculation gives G s « 0.1559 Tj J and 
G z « -0.1686 Tj J, which after substitution in Eq. © 
leads to SF » 7T 2 / J, where 7 » 0.8 • 10~ 4 . Addition to 
Eq. JSJl of the terms associated with more distant pairs 
of bonds leads only to a slight reduction of the numerical 
coefficient to 7 w 0.7 ■ 10 -4 . 

Thus we have demonstrated that SF cx T 2 , as in 
the case of the antiferromagnetic XY-mode\ with a 
kagome lattice |l2l llflj . In situations, when so-called 
" order- from-disorder" mechanism \\A llfij works already 
at the harmonic level, the free energy difference be- 
tween the accidentally degenerate ground states of frus- 
trated XY models is proportional to the first power of 

The fluctuation induced free energy (per unit length) 
of a zero-energy domain wall on the background of a 
striped state (see Fig, lb) is also given by SF. It has 
been shown in Ref. |l9j that in the frustrated XF-models 
with accidental degeneracy which manifests itself in the 
possibility of formation of zero-energy domain walls, the 
temperature of the phase transition associated with vor- 
tex pattern disordering (that is with the proliferation of 
such walls) can be estimated as 

T c ^E K /]n[T c /SF(T c )}, (8) 

where Ek ~ J is the energy of a kink on a domain wall 
and 8F(T) is the fluctuation induced free energy of a 
domain wall (per unit length). For 7 ~ 10~ 4 the loga- 
rithmical factor in Eq. © is close to 12, which means 
that the extreme smallness of SF(T) in the considered 
problem leads to the reduction of T c (in comparison with 
Ek) by one order of magnitude. 

However, the extreme smallness of 5F(T) manifests it- 
self much more strongly in the relative prominence of the 
finite size effects. For SF(T) w jT 2 /J the (normalized) 
probability to have a domain wall crossing a finite sys- 
tem of the width L can be estimated as exp[— (7T/ J)L] 
and is much smaller than one only for L L c — J/jT, 
which in our case gives L c > 10°. 

This means that although anharmonic fluctuations in 
the fully frustrated XY model with a honeycomb lattice 
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lead to the existence of a thermodynamic phase transi- 
tion (related with vortex pattern ordering) at not too 
small temperature T c ~ 10 _1 J, the sizes of Josephson 
junction arrays available experimentally, as well as the 
sizes of the systems which can be simulated in numerical 
experiments are currently not sufficient for observation 
of this ordering. 

To conclude, in the present work we have demonstrated 
that the fully-frustrated XY model with a honeycomb 
lattice is the first example of a statistical model in which 
the accidental degeneracy of different ground states is 
not removed by the free energy of harmonic fluctuations 
(neither in the classical, nor in the quantum version of the 
model) , although in different states these fluctuations are 
described by different Hamiltonians. The responsibility 
for that can be traced to a hidden gauge symmetry, which 
manifests itself when these Hamiltonians are rewritten in 
terms of the mixed representation. 

We also have shown that accidental degeneracy is re- 
moved by the free energy of anharmonic fluctuations, 
which leads to the selection of a striped state of the 
type shown in Fig. la. The state with analogous vor- 
tex configuration is selected as well in a fully-frustrated 
superconducting wire network with honeycomb geometry 
(in the vicinity of the phase transition), but for entirely 
different reasons related with the possibility of the order 
parameter modulation [2l|. 

However, the estimates based on numerical calculation 
of the anharmonicity induced domain wall free energy 
show that the system size (L ^> 10 5 ), which is required 
for the observation of vortex ordering in a fully frustrated 
Josephson junction array with a honeycomb lattice is 
much larger than those which are typical for experiments 
or numerical simulations, which makes the observation of 
such an ordering rather problematic. It should be empha- 
sized that although the fully frustrated XY model with 
a honeycomb lattice has been the subject of Monte-Carlo 
simulations of Shih and Stroud |9j,|22j, these authors have 
not analyzed the structure of vortex pattern. 

The ideas developed here may also be extended to 
other geometries under current investigations. One of 
the most intriguing systems in this respect is the fully 
frustrated XY model on a dice lattice, which exhibits 
a similar degeneracy between its classical ground states 
[23|], and has been the subject of recent experiments [24| 
™d numerical simulations M- ^ particular . one of the 
main reasons for the absence of vortex ordering in mag- 
netic decoration experiments on Josephson junction ar- 
rays 13 , as well as in numerical simulations of Ref. [25] 
is very likely to be a not sufficent system size. 

The conclusion on relative prominence of finite size 
effects (leading to the destructiion of long range order) 
in situations when the stability of a vortex pattern is 
induced only by anharmonic fluctuations is applicable 
to even wider class of models. Their number includes 
the antifcrromagnetic XY model with a kagome lattice 



[I2I [Til Ejj. in which, in the thermodynamic limit, the 
anharmonic corrections to free energy lead to the stabi- 
lization of so-called V3 x V3 state at T < T c ~ 10 -4 J 
[l3|. However, the estimates analogous to those per- 
formed above allow one to conclude that the observation 
of such an ordering is possible only when L ^> L c ~ 10 7 , 
that is only in truely macroscopic system. In numerical 
simulations of Ref. 26 this condition definitely was not 
satisfied. 
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